Fluctuation of the Correlation Dimension and the Inverse Participation Number at 

the Anderson Transition 
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The distribution of the correlation dimension in a power law band random matrix model having 
critical, i.e. multifractal, eigenstates is numerically investigated. It is shown that their probability 
distribution function has a fixed point as the system size is varied exactly at a value obtained from 
the scaling properties of the typical value of the inverse participation number. Therefore the state- 
to-state fluctuation of the correlation dimension is tightly linked to the scaling properties of the 
joint probability distribution of the eigenstates. 
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It is already well established [Q, |[ |) that the eigen- 
states at the localization-delocalization transition driven 
by disorder are multifractals characterized by a set of 
generalized dimensions, D q . This is reflected both in the 
statistical properties of each eigenstate and in the joint 
probability distribution of the eigenvector components. 
In the first case the multifractal analysis is performed on 
each eigenstate using the usual box-counting algorithm, 
while in the second case the scaling properties of the dis- 
tribution function are analyzed. 

The generalized dimensions are associated to the scal- 
ing of the qth moment of the eigenstates, ip a (r). These 
moments are also termed as generalized inverse partici- 
pation numbers Jj], ^ 
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where L is the linear dimension of the system. The 
inverse of / 2 , for example, roughly equals the number 
of nonzero wave function components, therefore it is a 
good and widely accepted measure of the extension of 
the states. For delocalized states J 2 ~ 1/N, while for 
localized ones I 2 ~ 0(1). At criticality it scales with L 
with exponent D2 (|l|). This scaling may be found for 
each state ipa (r) but then the corresponding exponents 
fluctuate. 

Another property of multifractality is related to the 
probability overlap of eigenstates with energy separation 
substantially exceeding the mean level separation. In 
particular D 2 describes these density correlations: hence 
the name of correlation dimension. It has been shown 
that for multifractal eigenstates these correlations decay 
slowly Q no matter how sparse these states are. The 
role of Z?2 is even more important as it has been related 
to the level compressibility at criticality [|[ || that is a 
remarkable relation that has been corroborated at least 
for the case of weak multifractality. 

The multifractal properties of the wave functions at 
criticality also shows up in the anomalous spreading of 
a wave packet, where the value of D2 plays an essential 



role §. 

The importance of the generalized dimensions and es- 
pecially of D2 therefore is obviously hard to overestimate. 
This is the reason why it came as a surprise when in 
H the authors demonstrated that the correlation dimen- 
sion, D 2 , has a relatively broad distribution even with 
increasing systems size of a three dimensional Anderson 
model. The authors concluded that this broadness re- 
mains in the thermodynamic limit and therefore there is 
no unique value of D2 in this limit. Such a statement 
would imply serious objections against one-parameter 
scaling theory which is corroborated by many other tech- 
niques . Based on numerical simulations on a random 
matrix model having eigenstates that exhibit multifrac- 
tality, however, it has been shown in j|, [n| that the fluc- 
tuations of I2 can indeed be characterized by a universal 
probability distribution function (PDF), II(/ 2 ) once I2 is 
normalized with the typical value, J 2 = -^/J^, where 
exp(lnJ 2 ): 
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n(in/ 2 ) = n JV (in h/iT)- 



(2) 



So the PDFs are simply horizontally shifted by 

(ln/ 2 ) oc D 2 \nN (3) 



ln/^ p 



to put on top of each other and obtain the universal curve 
of II(ln/ 2 ) showing that there is indeed a unique value of 
E>2 to characterize the scaling of / 2 . In the present work 
we wish to reconcile the apparent contradiction existing 
about the fluctuations of the correlation dimension D 2 
and to show how to correct the misinterpretation of Q 
and obtain results that are compatible with H [To| |. 

To this end we performed numerical simulations on 
the power law band random matrix ensemble that has 
been recently intensively investigated || [l0| O, p"2| . 
This model is more suitable for a quantitative analysis 
since it allows to obtain a wider range of linear extension 
[N = L) of the system as compared to the usual Ander- 
son model defined over an N = L 3 lattice of linear size L. 
Moreover, the model is capable to show many of the fea- 
tures of the traditional Anderson transition ^ |^, [l2| . We 
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FIG. 1: Distribution function of the inverse participation 
number n(ln/2) vs the rescaled value of \n(l2/I 2 yp ) for dif- 
ferent matrix sizes N. Open (filled) symbols correspond to 
(3 = 1 (fl = 2). The inset shows the scaling of the typical 
value, ln/^ p = <ha/ 2 ) vs ln(iV). 



FIG. 2: Distribution function of the correlation dimension 
V(T>2) for the random matrix model with different matrix 
sizes JV for the case of conserved time reversal symmetry, 
(3 = 1. 



will show that indeed the PDF of the correlation dimen- 
sion of the states is a broad function which eventually has 
a fixed point at a value very close to D 2 obtained from 
the scaling of I 2 yp '. Moreover, from the average value, 
(D 2 ), one can extract the same answer as well. 

The elements in these ensembles of N x N complex her- 
mitian matrices are Gaussian distributed random vari- 
ables with zero mean, (Hi j) = 0, and a variance 



which yields 



((H id Y)=[l + (\i-j\/b) 2 } *x 



1/2, 

V/3, 



i 7^ 3 
i = 3 



(4) 



where (3 = 1 or 2 depending on the global symmetry 
of the system. In the presence of time reversal symme- 
try, = 1, and the matrices are real symmetric, in its 
absence, (3 = 2. The parameter b is an effective band- 
width that serves as a continuous control parameter over 
a whole line of criticality || . In this work we set b = 1 as 
this proves to give the best convergence towards a scale 
invariant distribution function |J. However, at the end 
we will give the 6-dependence of the results as well. 

The states have been obtained using standard diago- 
nalization routines for the cases of N = 128, ... , 2048. 
The properties of the states in the vicinity of the middle 
of the band have been collected. We chose 25% of the 
states around E = 0, and verified that by narrowing this 
energy window the results do not alter but obviously the 
statistics becomes worse. 

Each individual eigenstate may show scaling with dif- 
ferent exponents D 2 . On the other hand the distribution 
of the IPN, I 2 ~ N D ' 2 , of the eigenstates shows univer- 
sal scaling properties ||, thus a relation between these 
properties is expected. The distribution function of D 2 
can be connected to that of the 1% (0) as 



V N (D 2 )dD 2 =U N (I 2 )dI 2 



(5) 



V N (D 2 ) = n N (N D2 )N D * ln(iV) 



(6) 



We expect the distribution of D 2 to tend to a univer- 
sal distribution as well as that of I 2 . Evers and Mirlin 
|| have proved that in particular I 2 is a self-averaging 
quantity and numerically confirmed that for the case of 
the present model. The same has been shown for the case 
of the Anderson model, too j| . Therefore as N — > 00 
we expect that the distribution 



V N {D 2 ) -» 5(D 2 - D 2 ). 



(7) 



This is in direct contradiction to the conclusions of m. 

Furthermore using Eq. (^) for V(D 2 ) and Eq. (|^) we 
can readily derive expressions of the mean value and the 
variance of D 2 in relation to the mean value and variance 
value of In I 2 as 



and 



. . ~ const 

(D 2 ) -D 2 = > 

V 1 In AT 



var(ln/ 2 ) 
var(L> 2 ) = — — tttt^ ► 0, 



(In iV) 5 



(8) 



(9) 



where the convergence to zero for the variance is ensured 

by @ 



var(ln/ 2 ) = cr^ - aN y , 



(10) 



with the irrelevant scaling y = D 2 /(2(3) > and also 
a > and a 00 > 0. These logarithmic convergences, Eqs. 
(||) and (H), seem to be too difficult to detect for the case 
of the traditional Anderson model [B, even though, Ref. 
0, showed the scaling @ very clearly. 
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FIG. 3: Distribution function of the correlation dimension 
V(D<z) for the random matrix model with different matrix 
sizes N for the case of broken time reversal symmetry, = 2. 

In Fig. |l| first we show that indeed the scaling of the 
PDF n^ln^) holds and obtain the value of D 2 from 
the typical value, (In 22) for both (3 = 1 and 2. Their 
values are listed in Table |. As a next step we turn to 
the PDF of the correlation dimension, D 2 . The corre- 
lation dimension of the individual eigenstates have been 
obtained using the standard box-counting algorithm jl). 
In Figs. H and || we show how P(D 2 ) evolves by doubling 
N . These figures give the same information as Fig. 3 of 
H for the case of the power law band random matrix en- 
semble. By increasing TV the distribution gets narrower, 
moreover, there is a fixed point right in the interval of 
0.7 < D 2 < 0.8 for (3 = 1 and at 0.85 < D 2 < 0.9 for 
(3 = 2. In the insets of both figures we plot this region 
enlarged in a semi-log plot. The fixed point clearly ap- 
pears at D* « 0.75 for /3 = 1 and D* » 0.87 for (3 = 2. 
The first value is precisely the one obtained via the scal- 
ing of the typical value (ln/ 2 ) in ||] and is close to our 
numerical results in Fig. [I]. 

We argue that due to the self-averaging property of I 2 
H V(D 2 ) should tend to a ((-function so for finite N, N' 
there should be a fixed point in the distribution function, 

V N (D*)=V N ,(D*), (11) 

provided D 2 falls within the set of D 2 s of the eigenstates. 
This ensures that D 2 = D 2 is the unique value of the cor- 
relation dimension at N — > 00. Moreover it is also the 
average value, (D 2 ) in this limit (ph. The convergence of 
(D 2 ) towards D 2 is given in Fig.^J. The value of (D 2 ) 
extracted from this fit together with the other estimates 
of D 2 are listed in Table Q. The error bars characterize 
the goodness of the linear fits as depicted in Figs. |l] and 
|]. The different approaches give slightly different values, 
some of them coincide within a few percent difference. 
We attribute the discrepancies to possible finite size ef- 



FIG. 4: Scaling of the mean value of the correlation dimen- 
sion as (D2) ln(iV) vs ln(iV). The lines are linear fits to the 
data. According to (^) the slop corresponds to D2 describing 
the scaling of \nll yp . 



TABLE I: Correlation dimensions obtained using different 
approaches 





= 1 


(3 = 2 


D2 


0.7036 ± 0.0024 


0.8576 ±0.0018 


(D2) 


0.7598 ± 0.0005 


0.9046 ± 0.0011 


D* 2 


0.7519 ± 0.0004 


0.8781 ± 0.0014 



fects and to the fact that the number of samples has not 
been large enough. For conserved time reversal symme- 
try, (3 = 1, we observe a coincidence between our (D 2 ) 
and D 2 values, together with the one obtained by M over 
a larger ensemble of samples and a wider range of ma- 
trix size. All the values from different estimates seem to 
be within a range of 1-8% difference. Therefore we may 
conclude that despite the relatively large state-to-state 
fluctuations of the correlation dimension depicted in a 
wide PDF of V(D 2 ) in Figs. | and | one can still define 
a unique value of D 2 that is characteristic to the critical 
point. 

Let us apply a similar argument to the results pre- 
sented in Ref. ]||. Inspecting Fig. 3 of that work 
closely, contrary to what is claimed there, the distribu- 
tion does not show a universal, broad form but rather 
gets more and more peaked as N increases. Even the 
authors admit a very large drop of V(D 2 ) in the interval 
of 1.5 < D 2 < 2.0 that apparently is system size depen- 
dent. One may find a similar fixed point as in our Figures 
|l| and H at around D 2 « 1.55. This value is consistent 
with most of the presently known estimates of the corre- 
lation dimension calculated using different methods and 
listed in It is, however, somewhat larger then more 
recent and accurate estimates J3|, [l0|. The possible rea- 
son for this discrepancy has been discussed in detail in 
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FIG. 5: The b/3 dependence of the correlation dimension D2 
obtained from its average value. The dashed lines correspond 
to the analytical estimates in the limiting cases of 6 <C 1 and 
6» 1 ©. 



Finally we present the behavior of the correlation di- 
mension as a function of the parameter b. This is pre- 
sented in Fig. (||) which shows that in fact D 2 is a unique 
function of bp. In this figure we plotted the slope of 
(D2) In A vs IniV for the various values of b and (3. The 
analytical estimates of || 



D 2 = 



1 - (7T&/3)- 

26/3, 



6» 1 
6< 1. 



(12) 



According to the figure we see that our numerical simula- 
tions are in full accordance with the analytical estimates. 
The discrepancies here also are attributed to the consid- 
erably smaller set of samples and smaller range of matrix 
size N as compared to Q. 

In conclusion we have performed numerical simulations 
of the power law band random matrices in order to ob- 
tain the PDF of the correlation dimension. We found 
that this function has a fixed point, D^, very close to 
the value of D2 obtained from the scaling of the typical 
IPN, ll vp . Besides a direct relation between the average 
correlation dimension, (D2), enabled us to obtain a third 
estimate of the correlation dimension. All of these are 



in good agreement with each other that corroborates the 
expectation, that there is a unique value of the correla- 
tion dimension characterizing the states at the Anderson 
transition. For finite size simulations, the state-to-state 
fluctuations of both the IPN and the correlation dimen- 
sion are large, nevertheless, both the PDF of I2 and that 
of T>2 provide a unique value of D2 in the limit N — > oo. 
This way we have reconciled the apparent contradiction 
between the state-to-state fluctuation of the correlation 
dimension and the universal scaling of the IPN. Similar 
results may be achieved for other generalized dimensions, 
e.g. the information dimension, D\ 
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